Torsion Cosmology and the Accelerating Universe 
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Investigations of the dynamic modes of the Poincare gauge theory of gravity found only two good 
propagating torsion modes; they are effectively a scalar and a pseudoscalar. Cosmology affords a 
natural situation where one might see observational effects of these modes. Here we consider only 
the "scalar torsion" mode. This mode has certain distinctive and interesting qualities. In particular 
this type of torsion does not interact directly with any known matter and it allows a critical non-zero 
value for the afiine scalar curvature. Via numerical evolution of the coupled nonlinear equations 
we show that this mode can contribute an oscillating aspect to the expansion rate of the Universe. 
From the examination of specific cases of the parameters and initial conditions we show that for 
suitable ranges of the parameters the dynamic "scalar torsion" model can display features similar 
to those of the presently observed accelerating universe. 

PACS numbers: 98.80.-k, 98.80.Jk, 04.50.+h, 98.80.Es 



I. INTRODUCTION 



One of the outstanding successes of theoretical physics 
in the latter part of the last century which led to a much 
deepened understanding was the recognition that all the 
known fundamental physical interactions, the strong, 
weak, and electromagnetic — not excepting gravity — can 
be well described in terms of a single unifying princi- 
ple: that of local gauge theory. Although there are other 
possible gauge approaches, for gravity it seems highly ap- 
propriate to regard it a gauge theory for the local symme- 
trygroup of Minkowski space time: the Poincare group 
[l|, Q . Such a consideration led to the development of the 
Poincare Gauge Theory of gravity (PGT) [l,Ei,i,0,i- 
The PGT has a priori independent local rotation and 
translation potentials, which correspond to the metric- 
compatible connection and orthonormal co-frame; their 
associated field strengths are the curvature and torsion. 
The spacetime then has generically a Riemann-Cartan 
geometry. Because of its gauge structure and geometric 
properties the PGT has been regarded as an attractive 
alternative to general relativity. The general theory in- 
cludes as exceptional cases Einstein's general relativity 
(GR) with vanishing torsion, the Einstein-Cartan theory 
with non-dynamic torsion algebraically coupled to the in- 
trinsic spin of the source, as well as the teleparallel the- 
ories wherein curvature vanishes and torsion represents 
the gravitational force (a sort of opposite to Riemannian 
geometry). Aside from these exceptions the generic PGT 
has, in addition to the metric familiar from Einstein's 
GR, a connection with some independent dynamics. This 
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additional dynamics is reflected in the torsion tensor. 

There is a natural physical source for the torsion of 
spacetime: namely spin 1/2 fermions. The effect is gener- 
ally assumed to be small at ordinary densities, but could 
have a major influence at high densities (e.g., beyond 
10'^* gm/cm'^), and thus it was expected to have impor- 
tant physical effects in the early universe Torsion 
cosmology investigations were initiated by Kopchyski . 
Some early investigations attracted attention especially 
because they noted that torsion might prevent the (at 
that time newly recognized) singularities. However, this 
hope quickly faded. Indeed it soon was argued that non- 
linear torsion effects were more likely to produce stronger 
singularities 

The various PGT dynamic modes beyond those of 
the metric were first investigated via the linearized the- 
ory (for outstanding examples of such investigations see 
0, El)- To this order the connection dynamics (which 
can be represented by the torsion tensor) decomposes 
into six modes with certain spins and parity: 2^, 1^,0*. 
Many possible combinations of well behaved (carrying 
positive energy at speed < c, criterion often referred to 
as "no ghost, no tachyon") propagating modes in the lin- 
ear PGT theory were identified. They were classified into 
about a dozen separate cases, almost any combination of 
up to 3 dynamic modes is allowed. Some nice investiga- 
tions of the PGT theory were also made using the Hamil- 
tonian analysis [1, [l^, i with findings consistent with 
the conclusions of the linearized investigation. Later, 
however, some potential problems were identified [14] . 
This prompted deeper investigations, which noted that 
effects due to non-linearities in the constraints could be 
expected to render most of the aforementioned dynamic 
cases physically unacceptable J^] . A fundamental inves- 
tigation identified two special cases, the so-called "scalar 
torsion" modes, which could be proved to be problem 
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free, having a well posed initial value prob lem [l6l | . Sub- 
sequently Hamiltonian investigations [l7l . [isj supported 
the conclusion that these two dynamic "scalar torsion" 
modes may well be the only physically acceptable dy- 
namic PGT torsion modes. 

In one mode (referred to as the "pseudoscalar" because 
of its 0^ spin content) only the axial vector torsion is dy- 
namic. (As a consequence of the dynamic field equations 
it turns out to be dual to the gradient of a scalar field; 
however it is not possible to treat this scalar field as the 
primary dynamical object without changing the nature of 
the theory [l^). Axial torsion is naturally driven by the 
intrinsic spin of fundamental fermions; in turn it natu- 
rally interacts with such sources. Thus for this mode one 
has some observational constraints [13] ■ Note that ex- 
cept in the early universe one does not expect large spin 
densities. Consequently it is generally thought that ax- 
ial torsion must be small and have small effects at the 
present time. This is one reason why we do not focus on 
this mode here. 

The other good mode, O"*", the so-called "scalar tor- 
sion" mode, has a certain type of dynamic vector tor- 
sion. (As a consequence of the dynamic equations given 
below in §11. B, it too turns out to be the gradient of 
a scalar field; this scalar field cannot, however, be re- 
garded as a fundamental potential — for essentially the 
same reasons as those mentioned in connection with the 
"pseudoscalar" mode). There is no known fundamental 
source which directly excites this mode. Conversely this 
type of torsion does not interact in any direct obvious 
fashion with any familiar type of matter [2ll |. Hence we 
do not have much in the way of constraints as to its mag- 
nitude. We could imagine it as having significant mag- 
nitude and yet not being dramatically noticed — except 
indirectly through the non-linear equations. This mode 
in particular has also attracted our interest because of a 
conspicuous consequence of the non-linear equations: in 
this case there is a critical non-zero value for the affine 
scalar curvature. 

Our theoretical PGT analysis thus led us to consider 
just two dynamic torsion modes. An obvious place where 
we might see some physical evidence for these modes is 
in cosmological models. The homogeneous and isotropic 
assumptions of cosmology greatly restrict the possible 
types of non- vanishing fields. Curiously, for torsion there 
are only two possibilities: 0"*", i.e., vector torsion which, 
moreover, has only a time component (and is thus effec- 
tively the gradient of a time-dependent function), and 
axial torsion, 0~, which is effectively the dual of a vector 
with only a time component (and thus can be specified 
as the gradient of a time-dependent function). Hence 
the homogeneous and isotropic cosmologies are naturally 
very suitable for the exploration of the physics of the 
dynamic PGT "scalar modes" . 

Thus cosmological models offer a situation where dy- 
namic torsion may lead to observable effects. Here we will 
not focus on the early universe, where one could surely 
expect large effects (although their signature would have 



to be disentangled from other large effects), and instead 
inquire whether one can see traces of torsion effects today. 
In particular we will here consider accounting for the out- 
standing present day mystery: the accelerated universe, 
in terms of an alternate gravity theory with an additional 
natural dynamic geometric quantity: torsion p^ . 

The observed accelerating expansion of the Universe 
suggested the existence of a kind of dark energy with a 
negative pressure. The idea of a dark energy is one of the 
greatest challenges for our current understanding of fun- 
damental physics [23, 24, 25]. Among a number of possi- 
bilities to describe this dark energy component, the sim- 
plest may well be by means of a cosmological constant A. 
However, there are some reasons for dissatisfaction with 
this model. In particular the so-called the cosmological 
constant problem notes that the theoretically estimated 
value of the vacuum energy density is about 10^^" times 
larger than the inferred cosmological constant. Moreover 
the coincidence or fine-tuning problem notes that it is 
highly unlikely that we should be living in the relatively 
short era when the rapidly changing ratio of the material 
energy and the cosmological constant is nearly unity. 

In the light of these problems there have been many 
interesting dynamical dark energy proposals. A popular 
idea is some unusual type of minimally coupled scalar 
field $ (quintessence field) which has not yet reached 
its ground state and whose current dynamics is basically 
determined by its potential energy V{^). This idea has 
received much attention over the past few years and a 
considerable effort has been made in understanding the 
role of quintessence fields on the dynamics of the Uni- 
verse (see, e.g., [26', [13, H^). However, without a specific 
motivation from fundamental physics for the light scalar 
fields, these quintessence models can be constructed rel- 
atively arbitrarily. There is a lot of room for speculation. 

Here we consider another possibility for explaining the 
accelerating universe: dynamic scalar torsion. We ex- 
plore the possibility that the dynamic PGT connection, 
refiected in dynamic PGT torsion, provides the acceler- 
ating force in the universe. As noted above, there are cer- 
tain "scalar torsion" modes which could have dynamical 
behavior. They could naturally provide the accelerating 
force in the universe. Here we will show that the effect 
of torsion can not only make the expansion rate oscillate, 
but also can force the universe to naturally have an ac- 
celerating expansion in some periods and a decelerating 
expansion at other times. Scalar torsion cosmology can 
avoid some of the problems which occur in other models. 

A comprehensive survey of the PGT cosmological mod- 
els was presented some time ago by Goenner and Miiller- 
Hoissen [1^. Although that work only solved in detail 
a few particular cases, it developed the equations for all 
the PGT cases — including those for the particular model 
we consider here. However that work was done prior to 
the discovery of the accelerating universe, and torsion 
was thus imagined as playing a big role only at high den- 
sities in the early universe. More recently investigators 
have begun to consider torsion as a possible cause of the 
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) but the subject 



accelerating universe (see, e.g., [30|, 
has not yet been explored in detail 

We have taken a first step in the exploration of the 
possible evolution of the Universe with the scalar torsion 
mode of the PGT. The main motivation is two-fold: (1) 
to have a better understanding of the PGT, in particu- 
lar the possible physics of the dynamic "scalar torsion" 
modes; (2) to consider the prospects of accounting for 
the outstanding present day mystery — the accelerating 
universe — in terms of an alternate gravity theory, more 
particularly in terms of the PGT dynamic torsion. With 
the usual assumptions of isotropy and homogeneity in 
cosmology, we find that, under the model, the Universe 
will oscillate with generic choices of the parameters. The 
torsion field in the model plays the role of the impercep- 
tible "dark energy" . With a certain range of parame- 
ter choices, it can account for the current status of the 
Universe, i.e., an accelerating expanding universe with a 
value of the Hubble constant which is approximately the 
present one. These promising results should encourage 
further investigations of this model, with a detailed com- 
parison of its predictions with the observational data. 

The remainder of this work is organized as follows: We 
summarize the formulation of the PGT and the "scalar 
torsion" mode in Sec. II, and then translate the equations 
into a certain effective Riemannian form in Section III; 
the specialization of these relations to the form describ- 
ing a cosmological model are presented in Sec. IV. Then a 
preliminary analytical analysis aimed at revealing the be- 
havior of the solutions is presented in Sec. V. In Section 
VI we present the results of our numerical demonstra- 
tions for various choices of the parameters and the initial 
data. The implications of our findings are discussed in 
Section VII and Sec. VIII is a conclusion. 
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0. 
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From the frame one constructs some auxiliary quantities: 
the reciprocal frame ^, which satisfies e'^e^'' = 5^'^ and 
eVcj^ = <^/i and the metric — ei^^ej^rj^j.^. Here our 
conventions are as follows: the Greek indices are the local 
Lorentz indices; whereas the Latin indices arc the coor- 
dinate indices. We use the metric signature (—,+,+,+). 

Following the standard paradigm, the conventional 
form of the PGT action, which is invariant under local 
Poincare gauge transformations, is taken to have the form 



A = 



j d*xe{LQ 



+ Lm). 



(5) 



Here e = det(e,^), eLG(e,'^, 9^6,'^, ^jF,^'') = 
eLQ{ei^ ,Tij^ , Rij'^") is the geometric gravity Lagrangian 
density and eLuiejT^ip, dip) — eLuiei^ Diip) is the 
minimally coupled source Lagrangian density, where ip 
represents all the matter and other interaction fields. We 
will not explicitly need the field equations for the non- 
geometric fields. Varying with respect to the geometric- 
gauge potentials gives the gravitational field equations. 
As explained in detail in the aforementioned references, 
they take the form 



E„ 



T * 



X7 ff — R * — ' 

with the field momenta 



deLc 



deLc 



(6) 
(7) 

(8) 
(9) 



and 



II. THE FIELD EQUATIONS 

A. Poincare gauge theory of gravitation 

Our considerations in this work are entirely classical. 
The form of the gravity theory we wish to consider here, 
the PGT, was worked out some time ago on the basis 
of the fundamental principles of gauge theory and ge- 
ometry [H, i, i, i, i, i, 0, 1. In the PGT there are 
two sets of local gauge potentials, the orthonormal frame 
field (tetrad) e^'' and the metric-compatible connection 
Tj/i'^j which are associated with the translation and the 
Lorentz subgroups of the Poincare gauge group, respec- 
tively. The field strengths associated with the frame and 
connection are the torsion 



T,/-2(a[,e,f'+F[,|.^e|,f), 
and the curvature 



(1) 



Rijf^" - 2{d[irj]^'' + Ffii^Tij]^'^), (2) 
which satisfy the Bianchi identities 



R 



(3) 



Ef,' := e^eLG-T^j'H,'' -R^.rH,a'\ (10) 
E^^' := Hi^^-\-. (11) 



The source terms here 



del 



M 



deLu 



(12) 
(13) 



are, respectively, the Noether energy-momentum and 
spin density currents, which (as a consequence of the as- 
sumed minimal coupling) automatically satisfy suitable 
energy-momentum and angular momentum conservation 
laws. 

The Lagrangian is chosen (as usual in gauge theories) 
to be at most of quadratic order in the field strengths, 
then the field momenta are linear in the field strengths: 



fc=i 



(14) 



k=l 
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here the three 



(fc).. 



(fc).. 

^ and the six R-'^ are the alge- 
braicaUy irreducible parts of the torsion and the curva- 
ture, respectively. The torsion in particular splits into 
the algebraically irreducible torsion vector, axial vector 
and tensor: 



density for this spin-O"*" mode is then 



ao "2 



ai 



AT^T^) , (18) 



T — T J 



where T„ 



The Hamiltonian analysis showed 



(16) 



which recompose to give 

4 2 1 „ 

Tijk — -^Q[i3]k + -^T[igj^k + -eijkmP" 



(17) 



The ak and 6fc in the Lagrangian are free coupling pa- 
rameters. Due to the Bach-Lanczos identity only five 
of the six 6^.'s are independent, oq is the coupling pa- 
rameter of the scalar curvature R := R^i,^^. Note that 
(because of the assumed quadratic Lagrangian, linear-in- 
field strength canonical momenta) one obtains, as in the 
standard physics paradigm, 2nd order equations for the 
potentials by varying the Lagrangian Lg of the PGT in- 
dependently with respect to the frame and connection. 
It should be remarked that these PGT equations are 
quite different from the problematical 4th order type of 
equations obtained from Riemannian geometry based La- 
grangians of the form R + (R,/ )^ when varied with re- 
spect to the metric. 

In the PGT, in addition to the dynamic metric repre- 
sented by the translational gauge potential (the orthonor- 
mal frame), the rotational gauge potential (the connec- 
tion) has some independent dynamics. As in other gauge 
theories, it is usually convenient to describe the dynamics 
of the connection (a non-covariant, gauge dependent po- 
tential) in terms of a tensorial field strength. In the PGT 
case these modes can be described by the torsion tensor. 
As mentioned in the introduction, the various PGT dy- 
namic torsion modes were first investigated via the lin- 
earized theory; it was shown that the torsion decomposes 
into six modes with certain s pins and parity: 2^, 1^,0^. 
Later investigations [H, [l6l. llTl [Tsj concluded that ef- 
fects due to non-linearities in the constraints could be 
expected to render all of these cases physically unaccept- 
able except for the two "scalar torsion" modes: spin-0+ 
and spin-0~. These two dynamic scalar torsion modes 
apparently are the only physically acceptable dynamic 
PGT torsion modes. 



B. simple spin-0^ mode 



that the number of degrees of freedom in Lg is three: the 
scalar torsion mode and two helicity states of the usual 
massless graviton (provided the scalar torsion mode is 
massive, i.e, ai ^ ao). It is necessary to impose certain 
sig n conditions on the parameters (see [j, M, llll, fisi - 



ai > , 



6 > 0. 



(19) 



There is a simple argument which accounts for the signs 
of these two parameters. In order to have least action 
the kinetic energy contribution from any dynamic vari- 
able must be positive (for if such a term were negative the 
action would have no lower bound, since we could have an 
arbitrarily large time rate of change for a dynamic vari- 
able). Consider that b is the parameter associated with 
the quadratic scalar curvature term R^. With the help of 
Eq. ([2]), it can be seen that the scalar curvature includes 
some time derivatives of one of our basic dynamic fields, 
the connection components: R^ — {e'^^e^ 



(2e* 



4(e* 



> 0. Hence 



the coefficient of this term in the action should be pos- 
itive. A similar argument based on ([T]), taking into ac- 
count the chosen metric signature and the restricted form 
of the torsion in our model as a consequence of the field 
equations, Eq. below, gives the sign of ai. 

Varying Lg (fT8| with respect to the potentials e^'^, r^'"^ 
gives gives the specific second order field equations of the 
general form (|6|7p for this mode. Assuming S^u^ = (i.e., 
the source spin current is negligible) we investigate first 
Eq. obtaining for this mode the three decomposed 
equations: 



= 

= 



(20) 
(21) 
(22) 



where /i := oi — oq is the effective mass of the linearized 
0^ mode. There is a special case with no dynamical 
scalar torsion if i? = — 6/Lt/6, a constant. We do not treat 
this exceptional degenerate situation as an isolated case 
(it is considered below as a limit of the generic case). 
Assuming that R + 6fi/b ^ generically leads to 



= 0. 



(23) 



Here we only investigate the simple spin-0^ case, i.e., 



choosing 02 — — 2ai, 03 
to vanish except for &(, = 
of this case please see [17 



= — ai/2 and taking all the &fc's 
- b 0. (For a detailed analysis 
.) Our gravitational Lagrangian 



Using these two constraints gives the restricted form of 
the torsion: 



T- Ai 

J- 1 J 



(24) 
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Substituting into Eq. ^ with our specific parameter 
choices gives the restricted field equation 

ML 2 ^ 24 3 ^ 

+R^'[^R-ao]}^T^\ (25) 

Now we have a complete set of the field equations, ((20)) 
and (l25|) along with ([23]) . In it was argued that this 
system had a well posed initial value problem. However 
that may be, the two main field equations are rather 
complicated. They really look nothing like the familiar, 
well-analyzed equations of GR. To help understand the 
significance of these new relations, and to use our previ- 
ous experience, we will do a translation of (|20l25p into 
a certain effective Riemannian form — transcribing from 
quantities expressed in terms of the orthonormal tetrad 
and connection Ti^"^ into the ones expressed in terms 
of the metric gjk and torsion Ty'^. Then we can com- 
pare the result with the more familiar field equations in 
general relativity. 

C. Translation 

As is well-known, the PGT affine connection can be 
represented in the form 

Ty ' = r,/ + i(T,/ + + T^-O , (26) 

where Tij'^ is the Levi-Civita connection, 

— i.9'tnj,i ^" 9rrti,j Qij.m) : (^'^) 

and Tij^ is the torsion. Accordingly the affine Ricci cur- 
vature and scalar curvature can be represented as 

R,^ = R,^ + V,T, + \{^k- Tfc)(T,,'= + T\, + T^O 

+ i(rfe™r'=™, + 2T,kmT^\) , (28) 

R = R + 2V^r + ^{T^jkT^'' + 2T^jkT'''' - ^T) , 

(29) 

where Rij and R are the Riemannian Ricci curvature 
and scalar curvature, respectively, and V is the covariant 
derivative with the connection Fy*^. 

For the case of interest here the torsion tensor has the 
restricted form (IM)) . Consequently the affine Ricci cur- 
vature and scalar curvature become 

R,j = i!„ + i(2Vjr, +5,,VfeT'=) 

+ ^m-g,,Tfcr'=), (30) 
R = R + 2ViT' - ^T,T' . (31) 



Applying this translation selectively [33*1 in Eqs. (|20|) and 
([25]) gives an alternate form of the field equations: 

V,i?+^(i?+^)r, -0, (32) 

aoiRij - -^gijR) - -R{R(ij) - -^gijR) 
-|^(V(,r_,-)-g,,VfeT'=) 
-g(2T,T, + g^.TkT'') - -T", , (33) 

while contracting Eq. with the help of Eq. (|HT|) yields 

aiR-nR = T. (34) 

Note that the relation (j33|l can be re-written into the 
form of Einstein's equation: 

ao{Rij - ^9ijR) = -Tij := -(TJj + Tlj) , (35) 

where Tlj is the source energy-momentum tensor and the 
contribution of the scalar torsion mode to the effective 
total energy-momentum tensor Tij is 

= - Y(V(,r,) - g., VfeT'^) - ^(2T,r, + g^.nr") 
- Ig^jR) ■ (36) 

However, it should be kept in mind that Tlj is only an 
effective quantity. Using this effective quantity allows 
us to use some of the insight we have obtained from our 
experience with GR. Thus we can regard the contribution 
of Tij to the rhs of ([35]) as something like that of an 
exotic field. This hybrid form is practical for our needs 
here, even though it is not really a proper fundamental 
physical description (one way to see this is to note that 
Tij cannot be obtained as the Hilbert energy-momentum 
density of some effective source Lagrangian). 

Eqs. p5l36p do allow us to appreciate some of the sim- 
ilarities and differences between this model and other ac- 
celerating universe models. However, to properly under- 
stand this model, one should consider the torsion dynam- 
ics geometrically rather than trying to regard it as just 
another field in a Riemannian spacetime obeying Ein- 
stein's equations. 

It is remarkable that the torsion vector — as a conse- 
quence of the field equation ([32]) — turns out to be the 
gradient of a scalar function. In fact we can identify 
the function as —(3/2) ln(i? + 6/i/6). However (especially 
given its geometric nature) we do not see any way to in- 
troduce this scalar potential directly into the Lagrangian 
as a fundamental field. Were that possible one could then 
directly compare features of our scalar torsion model with 
the various scalar field dark energy models. But as far as 
we can see, notwithstanding a few similarities, our model 
is really not much like those scalar field models. 
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III. FIELD EQUATIONS FOR TORSION 
COSMOLOGY 



For cosmology, assuming homogeneous and isotropic 
leads to the Friedmann-Robertson- Walker metric 



dr2 



1 — kr'^ 



) ,(37) 



For the effective energy-momentum tensor contribu- 
tion from the scalar torsion mode T^j, the explicit ex- 
pression is: 



r 



la — Id, 



1 



(51) 



where a{£) is the expansion factor, and fc is the curvature 
index. Here, to see the effects we are interested in as 
well as match the observations, it is sufhcient to consider 
only the simplest case: the flat universe with fc = 0. This 
yields the non-vanishing (Riemannian) Ricci and scalar 
curvature: 



and the off-diagonal terms vanish. 

We define pcft = P + Pt = — 3ao-ff^, where px = T^tt is 
the torsion-induced mass density, poff means the effective 
mass density which is deduced from general relativity. 
Pt = is an effective pressure due to contributions 
induced by the dynamic torsion. 



Rt = 3- 



3(H + H^ 



(38) 



R/ ^Rge^R/^^+2^ = H + , (39) 



a 



R = 6(- 



= 6{H + 2H^), 



(40) 



where H :— d/a. The torsion Ti should also be only time 
dependent, i.e., Ti=Ti{t). So from ([32|) the spatial parts 
of Ti vanish. Letting Tt{t)—^{t) we have 



Integrating this equation leads to 



R = 6{H + 2H^ - H<i>) - 2^ + -<^>^ 

3 



(41) 



.(42) 



From the field equations we can finally give the necessary 
equations to integrate: 



a 
H 

$ 
R 



aH , 

-1^R+J-T-2H\ 

Dfli OOi 

~-^R+7^T-3H'i> + l-<P\^ 
2ai 2ai 3 



(43) 
(44) 

(45) 

(46) 



where 
b 



(i? + ^)(3ff - <J>)2 - - 3a,H' ^Tu^p, 

(47) 
(48) 
(49) 



p = wp . 



IV. A PRELIMINARY ANALYSIS OF THE 
EQUATIONS 



Equations (|43| - [46)) are the main equations for the in- 
tegrations to evolve the system. Regarding the parame- 
ters in the field equations, the Newtonian limit requires 
Go = — (87rG)~^. We take ai > and 6 > to sat- 
isfy the energy positivity requirement . Moreover the 
no tachyon condition for the scalar torsion is then also 
satisfied: p, — ai — ao > (SttG)"^ > 0. 

Before the detailed results are shown, we briefly ana- 
lyze the equations to obtain some insight about their be- 
havior. Let us first study the behavior of the affine scalar 
curvature R. The second derivative of R with respect to 
time can be obtained by operating a time derivative on 
Eq. (gni) and using Eq. (|i5)) : 



R 



2 ■ 2, 6p- 



2aop 
aib 



R, 



(52) 



here we assumed that all the variables, i.e., H, R are 
much smaller than the coefficient of the leading order 
term, i.e., 2aQp/aib. For T, which appears on the rhs 
of (|^5)) . we know it consists at least of quadratic terms 
of H, R, and $ from (|T7)) . so it should be smaller than 
the other variables. This shows that the coefRcient of R 
on the right hand side of ([52]) is negative: 2aop/aib < 0. 
From this analysis we find that the late-time behavior of 
R will be essentially oscillating with the period 



T = 2tt, 



aib 
2aop 



(53) 



By a similar argument, it is easy to infer that $ has the 
same periodical behavior. 

Next we direct our attention to the behavior of the ex- 
pansion factor a. The acceleration of the expansion factor 
a can be obtained by combining Eqs. ()34|43l44|50l5ip : 



Here we consider only the matter-dominated era, where 
the pressure p is negligible. 



pR + T 

f 

6ai 



a 
a 



3pT + Pcff ^ 

6ao 



(54) 
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From this the relation between the acceleration of the 
expansion d and the quantity 3pT + PcS can be clearly 
seen. Since ao < 0, it shows that a > as long as 3pT + 
PcS < 0, and vice versa. We will discuss this relation and 
its demonstration in the next section. 

The period of a and H , if they exist, should be same as 
that of $ and R. Because the variables are all highly cou- 
pled to each other to form an equation set, there should 
generically exist a common period in the solution. This 
point will be demonstrated in the later numerical analy- 
sis. 

We need to look into the scaling features of this model 
before we can obtain the sort of results we seek on a 
cosmological scale. In terms of fundamental units we can 
scale the variables and the parameters as 



t 



t/i, a —t a, H — 



ao 



en, $ ^ £$, fR, 

-^fai, p.-^ffi, b^b, (55) 



where I = vSttG. So the variables and the scaled pa- 
rameters flo, ai, and b all become dimensionless, and 
ao = —1. Furthermore, Eqs. (|43H46p remain unchanged 
under such a scaling. However, as we are interested in 
the cosmological scale, it is practical to use another scal- 
ing to turn the numerical values of the scaled variables 
"gentler" (i.e., not stiff) from the numerical integration. 
In order to achieve this goal, let us introduce a dimen- 
sionless constant To, which represents the magnitude of 
the Hubble time (To = H^^ = 4.41504 x 10^^ seconds.) 
Then the scaling is 



t^Tot, a^a, H ^ H/Tq, $ $/To, R R/T^, 
ao -> ao, ai -> ai, ijl, b^ T^b, (56) 



With this scaling, all the field equations are kept un- 
changed while the period T TqT. 



A. parameter choice with constant scalar curvature 

Equation is of special interest among the field 
equations because of the existence of a constant value 
6/i/6. It shows that the scalar affine curvature remains 
a constant R = — 6/z/6 forever as long as its initial data 
has this special critical value. It is tempting to see how 
the system evolves with R = —6fj,/b initially. 

As mentioned in Sec. [Hi the positivity of the kinetic 
energy in the Hamiltonian analysis of the spin-O"*" case 
requires 6 > 0, oi > 0, thus ^ > since ao < 0. With 
such an assumption, the scalar affine curvature should 
not have the value R = —6fj,/b initially since this initial 
choice will require the matter density p to be negative 
from Eq. ([T7|) . Such a choice violates the assumption of 
energy positivity. 

However, if we tentatively relax the parameter require- 
ment for positive kinetic energy, i.e., allowing ai = —di < 
such that fi = —m < 0, the scenario will turn out to be 



quite intriguing. Under such a new parameter require- 
ment, if we set initially the scalar affine curvature 



R 



6p. 
T 



6m 



> 0, 



(57) 



then R will remain at this constant value for all the time. 
From Eqs. (|T7)l and ([50)1 . we can derive 



- 3aoH^ = p + pT > , 



where 



P 
Pt 



= 3diH^ 



2 6 



(58) 

(59) 
(60) 



Here the matter density will be positive as long as the pa- 
rameters are chosen suitably, such that diH'^ — vn? /2h > 
0. And the more interesting point is that the torsion- 
induced mass density px could "act like a dark energy" 
if the suitable parameter values are chosen. We can sim- 
plify the field equation ([44]) to 



H 



and it leads to 



o 2 

6 m 
Adib 



3 

2 M 



(61) 



(62) 



Combining with Eq. 



3TO'^/2ai6 > as long as 



3pT + Poff < 

show that 



0. By solving Eq. (pT|) . the solution will 



H 



V2II6 



for t oo. 



(63) 



By comparing this to a universe with a cosmological 
co nsta nt A where the Hubble function H approaches to 
■\/A/3 as i — > oo, we can see how to choose suitable val- 
ues such that the cosmological constant A and thus the 
dark energy can be mimicked in this torsion cosmological 
model with a constant afhne scalar curvature. We will 
demonstrate numerically the behavior of this case in the 
next section. 



V. NUMERICAL DEMONSTRATION 

In this section we would like to demonstrate two points: 
(1) The degenerate case R = —Qji/b with the relaxed pa- 
rameter choice mentioned in subsection IIV A| i.e., ai < 
and /i < instead of the normal choice ai > and p > 0. 
Although such a choice is against the positivity of kinetic 
energy, we would like to explore this scenario a little bit 
more, since it could mimic the cosmological constant and 
the other cosmological models with a negative kinetic 
energy [s^j- We can see that the torsion in the system 
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FIG. 1: Evolution of the Hubble function, H, the 2nd time derivative of the expansion factor, a, the temporal component 
of the torsion, and the afhne scalar curvature, R, as functions of time with the parameter choice and the initial data 
in Case I. The (black) solid lines represent the result of R{t = 0) = 6m/b, the (blue) dashed lines represent the result of 
R{t = 1) = 6m/b — 10~*, and the (red) dot-dashed lines represent the result of R{t = 1) = 6m/b + 10~* while all the other 
initial choices are fixed. 



becomes kinetic instead of being dynamic, and the ex- 
pansion is accelerating at late time; (2) In generic cases, 
i.e., R + 6fi/b ^ 0, with the proper parameter choice (i.e., 
ai > and /i > 0), the torsion in the system is dynamic, 
and its functional pattern has a periodic feature, i.e., it 
could be accelerating for a while, and then be followed 
by a period of deceleration with the pattern repeating. 
With suitable adjustments of the parameters and the ini- 
tial values of the fields involved, it is possible to change 
the period of the dynamic fields as well as their ampli- 
tudes. Furthermore, in the model, with some choices of 
the parameters and the initial values of the fields, it is 
possible to mimic the main apparent dynamic features 
of the Universe, i.e., the value of the Hubble function is 
the current Hubble constant in an accelerating universe 
after a period of time on the order of the Hubble time. 
In such a case, this model will describe an oscillating 



universe with a period on the order of magnitude of the 
Hubble time. This allows us to constrain the parameters 
and/or the value of the torsion field by comparing the 
observed data with the result from this model. 



The 4th-order Runge-Kutta method is applied for the 
integration of the field Eqs. The Universe is as- 

sumed to be matter-dominated, thus T « —p. The mass 
density p is determined from the fields via Eq. (|47|) . The 
fields and the parameters are scaled with Eq. ([55)1 and 
Eq. to be dimensionless, and to achieve a realistic 
cosmology. 
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FIG. 2; Evolution of the Hubble function, the 2nd time derivative of the expansion factor, the temporal component of the 
torsion, and the afhne scalar curvature as functions of time with the parameter choice and the initial data in Case II. 



A. Case I: constant R case 

In this case, the initial values of the fields are as follows: 

6m 



b ' 



a(io) = 50, H{to) = l, $(to) = 10, R{to) 
and the parameters are taken to be 

where to is the initial time: to = 1, the present time of 
our universe. Under this setting. 



~2b 



= 3 + 3m- 1.5 X lOW. 



(64) 



In order that the mass density in the current universe is 
about p « 30%, the parameter m is chosen as 

m = 0.012. 



This shows that H m/y/2dib w 0.84. The detailed 
result is shown in Fig. [1] where the evolved values of 
H, d, $, and R are plotted as the (black) solid curve in 
different panels. 

It is obvious in the bottom-right panel of Fig. [T] that 
the affine scalar curvature R remains constant, Qm/b. 
The behavior of the torsion $ can be understood through 
Eq. ([45]) . <& will increase (or decrease) until its value 
balances the rhs of Eq. (|15)) : this mainly depends on the 
sign change of the term 3-ff$ provided H > and T > 0. 
With the current initial choice in this case, decreases 
promptly at present until the balancing point is reached, 
as seen in the bottom- left panel of Fig.[TJ However, $ will 
not be a constant since the rhs of Eq. (j45|l still changes 
with time. The Hubble function H will always decrease, 
and approach to the fixed value m/ ^/2dib « 0.84, as 
shown in the upper-left panel of Fig. [H since the rhs 
of Eq. (pT|) is always negative. The acceleration of the 
expansion factor, a, is positive at late time, as seen in 
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the upper-right panel of Fig. [T] 

It is very interesting to see how the universe evolves 
if the scalar affine curvature R has a tiny deviation 
from the constant value Qm/h. Therefore, we chose 
the initial values of R as R{t = 1) = 6m/6 — 10~^ 
and R{t = 1) = Qm/h + 10^^ and evolved the system 
while keeping all the other initial choices the same as 
in the R = Gni/b case. The results are also plotted in 
Fig. [TJ In Fig. [U the (blue) dashed lines are for the 
R{t = 1) = Qm/b— 10~^ case and the (red) dot-dashed 
lines are for the R{t = 1) = 6m/b + 10~* case. The 
results show that once the scalar affine curvature R is 
smaller than Qm/b by a tiny amount, the values of R and 
thus the other fields will eventually return to a damped 
oscillating mode. Therefore, the universe will eventually 
approach a static condition. On the other hand, if the 
scalar affine curvature R is bigger than 6m/b by a tiny 
amount, the values of R and thus the other fields will rise 
unboundly. Either way the affine curvature will never re- 
cover its constant value. Therefore, the constant curva- 
ture case represents an unstable universe with an effective 
cosmological constant or a negative-kinetic-energy field. 
This phenomenon demonstrates the inherent instability 
of a system with a negative kinetic energy. However, as 
long as the deviation of R from 6m /b is small enough, 
it would be difficult from Fig. [T] to predict the future of 
the universe, since all the lines of these three cases are 
virtually overlapped together until a very late time. By 
careful fine tuning we can arrange for a large variety of 
outcomes. This "chaotic" behavior well illustrates just 
how we can lose all physical predictability if we allow 
such unphysical parameter choices. 

Although the above parameter choice has the virtue 
of explaining the accelerating expansion of the Universe 
and the cosmological constant, we cannot accept such a 
parameter choice here since it violates the fundamental 
assumption of the positivity of the kinetic energy. There- 
fore, in the following cases, we will return to our normal 
physical assumption, i.e., ai > and fi > 0. 

B. Case II: oscillating acceleration of a 

For this case, we take the initial values of the field to 

be 

a(0) = 10, H{0) 5 X IQ-^ 

$(0) = 2 X IQ--*, i?(0) = -2 X 10-3, 

and the parameters are taken to be 

H = l.2, 6 = 4, 

The results plotted in Fig. [5] show that a, $, and R are 
damped periodic. 

In particular R has a periodic character as shown in 
the bottom-right panel of Fig. [2] According to Eq. ([53]) 
its period is T = 2 tt -^—016/200/^ ~ 3.63, which is close 
to the period of the variables shown in Fig. [51 The most 



interesting part is the behavior of a, which is periodic 
with the same period as $ and R. As shown in the top- 
right panel of Fig. [21 d could be positive as well as being 
negative and the pattern of its function is similar to the 
pattern of R. Therefore the behavior of is a declining 
baseline plus a damped oscillation, as shown in the top- 
left panel of Fig. [3 

On a broader viewpoint of the evolution of this system, 
a, and R will be slowly damped, and H will approach 
zero after a long time. The important feature of this 
case is that the universe could oscillate due to dynamic 
torsion. In such a scenario the present day acceleration 
of the Universe is not so strange, d is oscillating and it 
happens to be increasing at this time. Furthermore, the 
oscillation period would be determined mainly by the pa- 
rameters and the initial values of the fields in this model. 
This encouraged us to try to find parameter values and 
initial conditions which more nearly resemble the current 
status of the Universe. Such a choice it will be shown in 
the next case. 



C. Case III: A presently accelerating universe 

In this case, we would like to compare the numerical 
values of the torsion model with the observational data 
of the universe. The initial data is set at the current time 
to — 1, after scaling, instead of to = 0. The parameters 
and initial conditions chosen are as follows 

a{to = 1) = 50, H{to = 1) = 1, $(to = 1) = 1.4, 

Y{to = 1) = i?(io = 1) + ^ = 6.2, 



and 

^i=l.09, 6=1.4. 

Here the initial data has been scaled according to 
Eqs. (|55l - [5S)) such that the current value of the Hubble 
function is unity. Therefore we get realistic values in our 
universe: the Hubble constant at present, H{to — 1), is 

H= ^ _.i«70J^. (65) 

4.41504 X 1017 s s • Mpc ^ ^ 

The results of the evolution with the parameters and ini- 
tial conditions are plotted in Fig. [3] In the top-left panel 
the Hubble function H is damped-oscillating at late time. 
In the top-right panel, it is obvious that d is damped and 
oscillating during the evolution and is positive at the cur- 
rent time t w 1, which means the expansion of the uni- 
verse is currently accelerating, ^{t) and R{t) are also 
plotted in Fig. [3] to show the correlation of the evolution 
between these variables. We observe that the values of 
the variables H{t), d{t), R{t), and $(t) become relatively 
high before t/Tg = 0.4. However this situation need not 
be taken too seriously, since it describes the earlier time 
of the universe, and our matter-dominated era assump- 
tion is not appropriate for such an early period of time. 
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FIG. 3: Evolution of the Hubble function, H , the 2nd time derivative of the expansion factor, a, the temporal component of 
the torsion, $, the afiine scalar curvature, R, the mass density, p, the effective mass density, p^s = p + pT, and the quantity 
3pT + PcS, as functions of time with the parameter choice and the initial data in Case III. 
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In order to have a deeper understanding of the settings 
of this case, the matter density p, the effective mass den- 
sity Peff = P + Pt, and the quantity 3pT + Pcff are plotted 
in the bottom panels of Fig. [31 The value of p, shown 
in the bottom-left panel, is decreasing at i « 1 as the 
universe is expanding and is always positive, while the 
effective mass density pofr, plotted in the same panel, 
shows an "oscillating" behavior around the curve of p. 
The oscillating behavior of peS comes from the contribu- 
tion of the torsion-induced mass density pT and simply 
indicates that px is not positive-definite in general. In 
fact the value of px turns from negative to positive when 
the time is around t « 0.7. As to the quantity 3pT + Peff, 
we can understand its importance for the evolution of 
the universe through Eq. (I54|) in which the value of this 
quantity decides the status of the acceleration. We can 
see this much more clearly by checking the correlation 
between the curves of a and 3pT + Peff in Fig. [31 Also by 
comparing the two bottom panels of Fig. [31 it is obvious 
that the torsion-induced pressure is negative when 
the universe accelerates, and positive when the universe 
decelerates. 

In this case the scaled value of p(t = 1) = 0.83 and 
its physical value is p{t = Tq) = 2.61 x 10~^°g/cm^. 
The Universe is supposed to be very close to the critical 
density, pc = 3c^H^/8ttG = 9.47 x lO'^Og/cm^; we find 
the ratio fi,„ = p/pc = 28%. In the standard ACDM 
model, flm ^ 30% with 5% baryonic matter and 25% 
dark matter. For our model fix = Pt/Pc = 72% acts 
like the energy density of dark energy. Therefore, this 
torsion model is able to describe a presently accelerat- 
ing expansion of the Universe with a proper amount of 
matter density. From the field equations we can see that 
the effect of the "dark energy" mainly comes from the 
nonlinearity of the field equation driven by the dynamic 
scalar torsion. 



Case 


p 


h 


H{1) $(1) 


y(i) 


a(l) 


a(i) 


P(l) 
10-«0g/cm^ 


III 


1.09 


1.4 


1 1.4 


6.2 


50 


27.59 


2.61 


IV 


1.27 


1.1 


1 0.8 


11.3 


50 


70.29 


5.23 


V 


1.38 


1.1 


1 1.1 


9.9 


50 


4.57 


2.48 



TABLE I: Here the parameter ao is set to be 1 in all of the 
three cases; -ff(l) means H{t = now), $(1) means <^{t — now), 
etc, under the scaling Eqs. (f55H56)l . 



supernovae data. Distance estimates from SN la light 
curves are derived from the luminosity distance 

'''-\/¥^"^"°<'><' +='/'!)• 

where Li^t and are the intrinsic luminosity and ob- 
served flux of the SN, and the rcdshift z = a{l)/a{t) — 1. 
Logarithmic measures of the flux (apparent magnitude, 
to) and luminosity (absolute magnitude, M) were used 
to derive the predicted distance modulus [35| 

p = TO-i\/ = 51ogiodL + 25, (67) 

where to is the flux (apparent magnitude), M is the lu- 
minosity (absolute magnitude), and c?l in the formula 
should be in units of megaparsecs. We found the rela- 
tions between the predicted distance modulus p and the 
redshift z in the three cases; they are plotted in Fig. [H 
For comparison, we also plot the prediction of the ACDM 
model with f2m = 0.3 and J^a = 0.7 by employing the fol- 
lowing formula js^ 



Jo v/(iTiHi~ra^^iF^i(2Ti)a^ ^ 

(68) 

The astronomical observational data [H, [s^l are also 
plotted in Fig. [H for comparison. The plots show that 
for small redshift z (e.g., z < 1.9) all three cases of the 
dynamical torsion models give an accelerating universe 
just like the ACDM model does. For larger z these cases 
might turn the Universe into a deceleration mode, which 
is consistent with the behavior of the various quantities 
shown in Fig. [31 We can see that Case V gives the clos- 
est curve behavior to the one from the ACDM model, al- 
though in Case V the matter density is only about 26% of 
the critical density. However, it was not meant to have a 
detailed comparison in this plot between our models with 
the ACDM model. Instead, in Fig.[4l we demonstrate the 
possibility of the scalar torsion field accounting for the ef- 
fect of dark energy with a suitable set of parameters and 
initial data. This allows us to study the dark energy 
problem from a new and different angle. 



VI. DISCUSSION 



D. Other Cases 

We continue to look at two more cases, which are listed 
in Table U along with Case III, obtained by taking differ- 
ent values of the parameters and the initial conditions, 
along with physical values of the significant mass density 
p. We find that the results of the other two cases have a 
behavior qualitatively similar to that of Case III. 

Now we would like to compare our results with the 



In this work we introduce into the evolution of a uni- 
verse without a cosmological constant a certain dynam- 
ical PGT scalar torsion mode taken from our earlier 
work p7j . From the assumption of the homogeneity and 
isotropy of the universe, only the temporal component 
of the torsion $ will survive and affect the evolution of 
the universe at late times. With the field equations (|i31 - 
I46p . we analyzed analytically and numerically the evo- 
lution of the system. We found that in generic cases, 
i.e., R + 6fj,/b ^ 0, with the proper parameter choice 
(i.e., fli > and p > 0), the torsion $ in the system 




FIG. 4: Comparison of different torsion models and tlie standard ACDM model with the observational data via the relation 
between the distance modulus /i and the redshift z. The supernovae data points, plotted with (brown) circles, come from [36| . 
the result of standard ACDM model (Sim = 0.3, JIa = 0.7) is plotted by the bold (red) dotted line. The results of Case III, IV, 
V, are represented by the bold solid line, the (blue) dashed line, and the (green) dot-dashed line, respectively. In the inset, the 
models and data are shown relative to an empty universe model (f2 = 0). 



is dynamic, and a, <i>, R tend to have a damped peri- 
odic behavior with the same period while the behavior 
of iJ is a declining baseline plus a damped oscillation. 
With certain choices of the parameters of /i and 6, and of 
the initial data of i?, <&, and i?, like Cases III-V in the 
previous section, this model can describe an oscillating 
universe with an accelerating expansion at the present 
time. 

Before we can give an adequate discussion of the via- 
bility of this model as an explanation of the accelerating 
universe, we should check whether this model can survive 
under the constraints of the theoretical and experimental 
tests. 

There have been numerous investigations on the ex- 
istence of torsion since this geom etric quantity entered 
the realm of gravity (see (sol. |40|. |4]| and the references 
therein). As mentioned above, this model has not only 



passed the important classical tests ("no-ghosts" and 
"no-tachyons" ) , it is also one of the two scalar torsion 
modes — the only PGT cases which are known to have 
a well posed initial value problem [iTj and which may 
well be the only viable dynamic PGT torsion modes that 
can evade the non-linear constraint problems. There 
have also been some laboratory tests in search of tor- 
sion [43 , |43| . The main idea among these experiments is 
the spin interaction between matter and torsion. The 
cosmological tests on torsion investigate the effect of 
torsion-induced spin flips of neutrinos in the early Uni- 
verse which could alter the helium abundance and have 
other effects on the early nucleosynthesis [13, . How- 
ever Dirac fermions interact only with the totally anti- 
symmetric pseudo-scalar part of the torsion. Thus these 
tests can only consider the pseudo-scalar mode (axial- 
vector torsion), not the scalar mode torsion used in our 



14 



model. The type of torsion used in our model does not in- 
teract directly with any known matter. Thus, these tests 
cannot really give a serious constraint on the amplitude 
of our scalar-mode torsion. 

Among the models in which torsion is applied to the 
cosmological problem, Capozziello et al. [3ll, have 
done a serious study on replacing the role of the cos- 
mological constant in the accelerating Universe. With a 
totally antisymmetric torsion without dynamical evolu- 
tion, their model is consistent with the observational data 
by tuning the amount of the torsion density, although 
this model cannot solve the coincidence problem. On the 
other hand, the oscillating universe models with a de- 
signed mechanism: an oscillating potential, an oscillating 
parameter of the equation of state, etc. [13, |H, aim 
to solve the coincidence problem. Here we found that our 
model takes some virtues from both kind of models, i.e., 
our model is capable of solving the coincidence problem 
of an accelerating universe with a dynamical scalar-mode 
torsion, which is naturally obtained from the geometry of 
the Riemann-Cartan spacetime, instead of from an exotic 
scalar field or a designed mechanism. 

If we consider the spacetime as Riemannian instead of 
Riemann-Cartan, by absorbing the contribution of the 
torsion of this model into the stress-energy tensor on the 
rhs of the Einstein equation, then this contribution will 
act as a source of the Riemannian metric, effectively like 
an exotic fluid with its mass density px and pressure 
varying with time (even though the time evolution of 
the torsion is not like that of a such a fluid). Moreover, 
the effective fluid appears to have presently a negative 
pressure, and consequently a negative parameter in the 
effective equation of state, i.e., wt = Pt/pt, which drives 
the universe into accelerating expansion. Note that there 
is no constraint on the value of o^t which appears here, 
and its value could vary from time to time. It should be 
stressed that this is not a real physical fluid situation; 
the truth is that wt is nothing like "a torsion field equa- 
tion of state" , it is just a proportionality factor between 
Pt and pt, two expressions which effectively summarize 
the contribution of torsion acting as a source of the met- 
ric. The ratio ujt is of interest only to help understand 
the acceleration of this model and to enable a limited 
comparison with other dark energy proposals. 

One might be concerned about the value of the param- 
eter b. Its value should be small enough to be consistent 
with the constraints on the affect of the quadratic order 
term on the large scale structure of universe. The 
values of b we choose, i.e., b/{aQTQ) in the conventional 
unit, are on the order of unity. These chosen values are 
bigger than the magnitude of a related parameter, esti- 
mated in [ssj; however one cannot expect that estimate 
to be applicable here — since in that work quadratic Rie- 
mannian curvature terms were considered (they lead to 
4th order field equations) instead of the afhne curvature 
terms we have used (which give 2nd order equations). 
As far as we know the parameter p does not have too 
much constraint on it, except for its positivity as a mass 



parameter, since the baryonic matter will only interact 
with the scalar torsion indirectly by gravitation. 

One may wonder: how large must the torsion be in or- 
der to produce observable effects in the the present day 
universe, e.g., the observed acceleration? conversely, how 
large can the torsion be without violating some observa- 
tional constraint? The questions merit a detailed study. 
Here is a simple argument that indicates a magnitude. 
Let us compare the terms in the Lagrangian density and 
the field equations for the PGT scalar torsion model and 
the Einstein theory with a cosmological constant. Note 
that the presumed cosmological constant is "so small" 
that it has no noticeable effect in the laboratory, nor on 
the solar system scale, nor on the galactic scale. Nev- 
ertheless it is large enough to have the dominant effect 
on the cosmological scale. Hence we are led to infer that 
we should consider aiT^ ~ ^ K ^ a^p ^ H'^ . With 
such a choice we can expect that torsion may be able 
to accelerate the universe and yet not be conspicuous on 
smaller scales. 

The 0+ torsion mode in this model effectively gives a 
scalar field, yet this scalar field is, in fact, quite different 
from the various scalar field models of "exotic matter" , 
e.g., the quintessence models, in several significant ways: 
(i) torsion cosmology is derived naturally from a geomet- 
ric gravitational theory, which is based on fundamental 
gauge principles, instead of on the hypothesis of the exis- 
tence of a dark energy tailored to producing an explana- 
tion of an accelerating universe; (ii) thus there are only a 
couple of free parameters in torsion cosmology, instead of 
an ad hoc potential that can be rather arbitrarily chosen 
to fit the observations. Therefore, a torsion cosmological 
model should be more restrictive, and should be easier 
to be confirmed or falsified; (iii) based on its tensorial 
character, the coupling of torsion to the other fields is 
nothing like that which has ever been advocated for hy- 
pothetical scalar fields. Consequently we see no way to 
simply replace the scalar mode torsion with an effectively 
equivalent quintessence model. Thus torsion cosmology 
and the quintessence models are characteristically differ- 
ent, even though there are some similarities. 

Due to its intriguing behavior, we also turned our at- 
tention to a degenerate case, R — —6p/b, in Case I with 
the relaxed parameter choice of ai < and ^ < in- 
stead of the normal choice ai > and p > 0. Although 
such a choice is against the positivity of kinetic energy, 
we explored this scenario since it could mimic the cosmo- 
logical constant and the other cosmological models with 
a negative kinetic energy. Indeed the result does show an 
accelerating universe at late time. However, our further 
numerical experiments also show that in such a case it 
describes a very unstable universe. A small perturbation 
from the constant curvature will cause a sudden change 
which would only become apparent at some time in the 
future. 
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VII. CONCLUSION 

In this work we considered the scalar torsion mode of 
the PGT on in a cosmological setting and proposed it as a 
viable model for explaining the current status of the Uni- 
verse. Besides having a better understanding of the PGT, 
we study the prospects of accounting for the outstand- 
ing present day mystery — the accelerating universe — in 
terms of an alternate gravity theory, more particularly 
in terms of the PGT dynamic torsion. With the usual 
assumptions of isotropy and homogeneity in cosmology, 
we find that, under the model, the Universe will oscillate 
with generic choices of the parameters. The torsion field 
in the model could play the role of dark energy. With 
a certain range of parameter choices, it can account for 
the current status of the Universe, i.e., an accelerating 
expanding universe with a value of the Hubble constant 
which is approximately the present one. Thus we have 
considered the possibility that a certain geometric field, 
dynamic scalar torsion — which is naturally expected from 
spacetime gauge theory — could fully account for the ac- 
celerated universe. 

The source of the torsion could come indirectly from 
the huge density of the particles with sufficient spin align- 
ment in the early universe. This scalar mode of torsion 
could be considered as a "phantom" field, at least in 
the matter-dominated epoch, since it will not interact 
directly with matter; it only interacts indirectly via grav- 
itation. Then the dynamics of the scalar torsion mode 
could drive the Universe in an oscillating fashion with 
an accelerating expansion at present. It is quite remark- 
able that a gauge theory of dynamic geometry naturally 
presents us with such a "phantom" field. This natural 
geometric field could act like a dark energy. 

However, there are also some points which need to be 
studied in much more detail before this model can more 
closely conform to reality. The model in Cases III-V of 
the previous section, suggests that the mass parameter of 
the torsion, /i, might be close to oq, and the parameter for 



the "kinetic" energy density of the torsion, b, may need 
to be as huge as Tq to achieve an accelerating universe. 
The restricted window of the parameter choices which 
allows a behavior like that of our universe might render 
the model less favored, even though the matter in the 
universe is not able to directly interact with the torsion. 
Meanwhile, the required choice of initial data and the 
values of the parameters may make this model unsuited 
to solving the fine-tuning problem. 

These dark sides should not be able to diminish the 
possibility of the scalar mode of the torsion in this model 
playing a significant role in the the evolution of the Uni- 
verse. The model has only a few adjustable parameters, 
so scalar torsion may be easily falsified — as "dark en- 
ergy." If it turns out that the accelerated universe can- 
not be explained in this way — that something else has the 
dominant dark energy role — it would still be reasonable 
to expect that there may be some observable cosmolog- 
ical effects from dynamic scalar torsion. Also, here we 
only used one of the viable modes of torsion in PGT; 
the our model will be more general if it is extended to 
include all the viable PGT torsion modes. We believe 
that future investigations along this line should be open 
to these possibilities. 
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